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We derive the effective low-energy theory for single-wall 
carbon nanotubes including the Coulomb interactions among 
electrons. The generic model found here consists of two spin- 
~ fermion chains which are coupled by the interaction. We 
analyze the theory using bosonization, renormalization-group 
techniques, and Majorana refermionization. Several exper- 
imentally relevant consequences of the breakdown of Fermi 
liquid theory observed here are discussed in detail, e.g., mag- 
netic instabilities, anomalous conductance laws, and impurity 
screening profiles. 



PACS numbers: 71.10.Pm, 71.20.Tx, 72.80.Rj 



I. INTRODUCTION 

The fascinating electronic and mechanical properties of 
carbon nanotubes have recently attracted a lot of atten- 
tion p]]. Nanotubes are tubular nanoscale objects which 
can be thought of as graphite sheets wrapped into a cylin- 
der. Shortly after their discovery during the carbon- arc 
fullerene production 0], theorists have made a surpris- 
ing prediction linking structural to electronic properties 
|^-^|. Specifically, the arrangement of carbon atoms 
on the tube surface is determined by the integer in- 
dices < m < n of the wrapping superlattice vector 
T = nai + too2, where a\ and ai are the primitive Bra- 
vais translation vectors of the honeycomb lattice. De- 
pending on the choice of n and m, a (n, m) nanotube 
should then be either a metal, a narrow-gap semiconduc- 
tor, or an insulator. This theoretical prediction has been 
amply confirmed in recent experiments [HJ . 

Nanotubes can be fabricated using either the carbon- 
arc method [p||r2[ or a novel laser ablation technique 
fiof , where Co- or Ni-doped graphite targets are laser- 
vaporized. The carbon-arc process usually yields multi- 
wall nanotubes (MWNT) composed of several concen- 
tric graphite sheets. One can then attach metallic leads 
to such a MWNT and measure, e.g., the magnetocon- 
ductance. The MWNT experiment by Langer et al. 
jl3| shows typical signs of a weakly disordered meso- 
scopic system, such as universal conductance fluctua- 
tions or weak localization phenomena. On the other 
hand, the laser ablation method allows for the controlled 
fabrication of single- wall nanotubes (SWNT). Rather 
large quantities of metallic (n, n) "armchair" nanotubes 
with n = 10 can be obtained using this technique [ {To) . 
In most cases, the SWNTs spontaneously crystallize in 



triangular-packed ropes containing « 20 to 100 SWNTs, 
to which one can again attach leads and perform trans- 
port measurements. The experimental results of Fischer 
et al. pd| have shown a linear temperature dependence 
of the resistance above a (non-universal) crossover tem- 
perature, with an increase of the resistance at lower tem- 
peratures. In addition, the experiments of Bockrath et 
al. [|l5| showed Coulomb charging effects 16 in such a 
rope. 

One can in fact also produce single SWNTs using the 
laser ablation process. From a fundamental (and the- 
oretical) point of view, a single SWNT is an intrigu- 
ing system, since Coulomb interactions induce a break- 
down of Fermi liquid theory [fl7| in any one-dimensional 
(ID) metal. Astonishingly, Tans et al. |l8| were able to 
attach contacts to a single 3/im long (10, 10) armchair 
SWNT. The transport measurements of Ref. |l^] were 
dominated by charging effects due to large contact resis- 
tances [around 500 kf2] between the leads and the SWNT. 
By finding a way to decrease the contact resistances one 
could circumvent charging effects, thereby allowing one 
to study the peculiar properties of a ID conductor. Nan- 
otubes are potentially much more stable ID conductors 
than conventional systems such as chain molecules [19|, 
edge states in the (fractional) quantum Hall effect [20|, 
or quantum wires in semiconductor heterostructures [21|. 
In these systems, non-Fermi liquid behavior has been no- 
toriously difficult to establish experimentally. 

As has been demonstrated in Ref. jig] , a metallic 
SWNT constitutes a perfect experimental realization of a 
ID conductor. Interacting ID electrons exhibit Luttingcr 
liquid rather than Fermi liquid behavior characterized by, 
e.g., the absence of Landau quasi-particles, spin-charge 
separation, suppression of the electron tunneling density 
of states, and interaction-dependent power laws for trans- 
port quantities. For detailed accounts of the physics of 
Luttinger liquids, see, e.g., Refs. 23|. In this paper 
we describe the nature of the non-Fermi-liquid state in a 
SWNT by deriving and analyzing the effective low-energy 
theory of carbon nanotubes. This state should be observ- 
able as soon as charging effects are overcome by lower 
contact resistances. A brief and incomplete account of 
our results has been given in Ref. J24J. 

The role of Coulomb interactions is most pronounced 
if the nanotube is metallic. What are the conditions 
for metallicity? A necessary condition arises because 
the Fermi vector should obey the transverse quantization 
condition Tk = 2ttI for an integer /. The first Brillouin 
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zone of the honeycomb lattice is a hexagon, and band- 
structure calculations |^5|,^6| show that the only gapless 
points are the corner points of this hexagon. One then 
finds that there are exactly two independent Fermi points 
denoted by K and K', with two linearly dispersing bands 
around each of these two points. Imposing the transverse 
quantization condition for the Fermi vector correspond- 
ing to K (or K') implies 2n + m = 31. If this condition is 
not fulfilled, the nanotube exhibits the "primary" band 
gap AE — 2v/3R « 1 eV, where v is the Fermi velocity 
and R the radius. We mention in passing that a nanotube 
without primary gap has a finite Fermi momentum only 
if n — m — 3dl, where d is the highest common divisor 
of n and m p7f . 

Even if the necessary condition 2n + m — 31 is fulfilled, 
the rearrangement of local bonds due to the curvature of 
the cylinder can introduce a "secondary" gap, AE » 10 
meV, which implies narrow-gap semiconducting behav- 
ior. For very small tube diameter, e.g., for a (6,0) tube, 
due to the strong curvature-induced hybridization of a 
and 7r orbitals, this effect can be quite pronounced |p8| . 
In the case of armchair nanotubes (n = m), however, the 
formation of a secondary gap is prevented by the high 
symmetry, and therefore a (n, n) SWNT stays metallic 
for all n. The highly symmetric structure of (n, n) [and 
(n, 0) "zig-zag"] nanotubes also prohibits chirality, i.e., 
the carbon atoms close around the waist of the nanotube. 
In contrast, all other nanotubes show chiral behavior [^9| . 

One-dimensional metals are expected to exhibit a 
Peierls instability due to the spontaneous formation of 
a lattice distortion. It has been argued in Refs. (^|| that 
this instability should not be of any practical importance 
in nanotubes. Mean-field estimates of the Peierls tran- 
sition temperature yield very low values (w 1 K), and 
fluctuations will then tend to wash out the instability 
even further. We therefore neglect coupling to lattice 
distortions in the following. 

So far we have discussed the case of a perfectly clean 
SWNT. Nevertheless, there are various sources for im- 
purities, e.g., structural imperfections of the nanotube 
(like substitutional atoms), charge defects in the sub- 
strate, topological defects, or twists. Topological defects 
are dislocations (kinks) that locally change the superlat- 
tice vector T by replacing one of the hexagons in the 
graphite network by a pentagon or a heptagon. Such a 
change typically introduces a metal-insulator interface. 
By combining a kink and its antikink (dislocation pair), 
it is possible to form local metal-insulator-metal junc- 
tions p0| . Transport is then dominated by tunnel events 
through this barrier. A less dramatic but probably more 
widespread source for backscattering is given by twists of 
the nanotube. These twists can originate from the mo- 
mentary position of the SWNT when deposited onto the 
leads, or from thermal fluctuation modes pl| . 

Generally, especially in the presence of impurities, pro- 
nounced effects of the Coulomb interaction on transport 



in ID metals can be expected [Q. As shown in Ref. f p3| , 
at very low energy scales transport is fully blocked by a 
single arbitrarily weak impurity. The conductance then 
vanishes as a function of temperature and/or voltage with 
anomalous interaction-dependent power laws. This ef- 
fect can be understood in terms of the Friedel oscillation 
p4] , |35| building up around a scatterer in a ID metal. 
The oscillatory charge screening cloud displays an alge- 
braically slow decay away from the impurity and thereby 
causes an important additional backscattering contribu- 
tion |36| . This crucial effect cannot be captured by linear 
screening p7| . The dielectric function approach to im- 
purity screening in SWNTs |38| is therefore practically 
useless in a determination of the Friedel oscillation. 

Below we focus on a single armchair SWNT, where 
interaction effects are very pronounced and direct exper- 
imental tests of the theory are within reach. However, 
as long as the above-mentioned band gaps are negligible, 
e.g., for sufficiently high temperatures, our effective low- 
energy Hamiltonian applies to any SWNT, even a chiral 
one. The only change arising for n m affects the pre- 
cise value of the interaction parameters introduced below. 
This remarkable circumstance is due to a U(l) symme- 
try of the graphite dispersion relation present close to 
the Fermi surface. Finally, we mention that by adding 
a bulk mass term, the properties of semiconducting or 
insulating SWNTs could be analyzed within the same 
framework p9| . 

The structure of this paper is as follows. In Sec. [] 
the effective low-energy description of SWNTs is derived. 
The resulting fermionic model is bosonized in Sec. |l|. A 
renormalization group analysis is carried out in Sec. IV 
and supplemented by the Majorana refermionization ap- 
proach described in Sec. [V| Experimentally relevant sus- 
ceptibilities and corre lation functions are discussed in 
Sec. In Sec. |VIl| , we study transport through the 
nanotu be, fo llowed by a discussion of impurity screening 

Ivni 



Sec. 



Finally, some concluding remarks can be 



found in Sec. IX. Throughout this paper, we put K = 1. 



II. EFFECTIVE LOW-ENERGY THEORY 

In this section we derive the effective low-energy de- 
scription of a single-wall (n, n) carbon nanotube. Our 
description applies at energy scales where the generic 
linear dispersion relation of a SWNT depicted in Fig. [j] 
is valid. For a (10, 10) nanotube, the theory therefore 
holds at room temperature and below. A similar ap- 
proach has been devised previously for the uncorrelated 
case by Kane and Mele |3l[] , building on earlier work |4C| ] 
for graphite intercalation compounds. Here we extend 
their theory and include the Coulomb interactions be- 
tween the electrons. Previously, this problem has only 
been investigated using the perturbative renormalization 
group (RG) for a weak short-range (Hubbard) interaction 
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pLE^. However, there is no external screening of the 
Coulomb interaction in the SWNT experiments of Ref. 
PH , and one has to take into account the long-ranged 
character of the Coulomb interaction potential (see also 
Rcf. Since one might employ suitable screening 

backgates in future experiments, the general (short- or 
long-ranged) case of an arbitrary interaction potential 
is treated here. The method used below allows us to 
get insight into the physics of the strong-coupling regime 
emerging at low temperatures. 



A. Low-energy expansion 

The remarkable electronic properties of carbon nan- 
otubes are due to the special bandstructure of the tt elec- 
trons in graphite. The simplest bandstructure calcula- 
tion starts from a nearest-neighbor tight-binding Hamil- 
tonian on the honeycomb lattice, which can be straight- 
forwardly diagonalized p5fl . Remarkably, the only gap- 
less points of the resulting dispersion relation are the cor- 
ner points of the hexagonal first Brillouin zone. Hence 
there are only two linearly independent Fermi points aK 
with a = ± instead of a continuous Fermi surface. Up to 
energy scales < 1 eV, the dispersion relation around the 
Fermi points is, to a very good approximation, linear, 



E±(q = k-K) = ±v\q\ 



(2.1) 



with the same relation around the other Fermi point. 
The + (— ) sign corresponds to the conduction (valence) 
band, respectively. Obviously, close to the Fermi surface 
a U(l) symmetry holds, since the direction of q does not 



matter in Eq. (2.1 



Since the basis of the honeycomb lattice contains two 
atoms, there are two sublattices p = ± shifted by the 
vector d — (0, d). with the nearest- neighbor C-C distance 
d = a/V3 = 1.42A (h ere a denotes the honeycomb lattice 
constant). As a consequence, there are two degenerate 
Bloch states 



(2.2) 



at each Fermi point a = ±, where r = (x,y) lives on 
the sublattice p under consideration. In Eq. (2.2), we 
have already anticipated the correct normalization for 
nanotubes. We follow Ref. |Q and choose the Bloch 
functions separately on each sublattice such that they 
vanish on the other. One can then expand the electron 
operator in terms of these Bloch waves (a = ± is the spin 
index) , 



pa 



(x,y) F pacr (x,y) 



(2.3) 



which introduces slowly varying operators F paa (x,y). 
The resulting second-quantized effective low-energy the- 
ory of a graphite sheet is given by the 2D massless Dirac 
Hamiltonian IHOfl, 



H G = 



P drF^ (tV) ,F , 

' j j ^' x paa\' v /pp'^p'. 



pp' aa 



(2.4) 



in accordance with the spectrum ( |2.l| ). Here r = (t x , t v ) 
are standard Pauli matrices. 

Wrapping the graphite sheet leads to the generic band- 
structure of a metallic SWNT shown in Fig. [l]. For a 
(n,n) armchair SWNT, the Fermi vector is K = (kp,0) 
with kp = 4ir/3a. We take the x-axis along the tube di- 
rection, and the circumferential variable is < y < 2irR. 
The armchair SWNT radius is R = \/3na/2w, which 
yields R — 1.38 nm for a (10, 10) nanotube. Quantiza- 
tion of transverse motion now allows for a contribution 
exp(imy / R) to the wavefunction. However, excitation of 
angular momentum states other than m = costs the 
energy w 10 eV/n. In an effective low-energy theory, we 
may thus omit all transport bands except m = and 
hence arrive at a ID situation with k y = 0. Since also 



y — 0, the Bloch states (2.2) do not depend on the 



circumferential coordinate, and the corresponding wave- 
functions are indeed stiff around the waist of the nan- 
otube. 

Instead of the low-energy expansion ( |2.3| ) , the electron 
operator is now written as |pl| 



pa 



ip pa (x 1 y)i> paa (x) , 



(2.5) 



which introduces ID fermion operators ipp a <r(x) depend- 
ing only on the x coordinate. They correspond to the 
Fp a a for the graphite sheet. Neglecting Coulomb inter- 
actions for the moment, the Hamiltonian can be read off 
from Fig. |l|, 



Ho = -v^p 



dx ip, 



(2.6) 



For this expression, we have chosen a preferred direction 
in Eq. ( |2.4| ). However, due to the mentioned U(l) in- 
variance, we can choose any direction without affecting 
the resulting low-energy Hamiltonian [we have explicitly 
checked that this holds even in the presence of interac- 
tions]. Switching from the sublattice (p = ±) description 
to the right- and left-movers (r = ±) indicated in Fig. |l| 
implies two copies (a = ±) of massless ID Dirac Hamil- 
tonians. By using a suitable gate, one can experimen- 
tally tune the average charge density on the nanotube 
and hence adjust the Fermi energy. In contrast to the 
half- filled band (Ep = 0) encountered in 2D graphite 
sheets, one is normally off half-filling in SWNT experi- 
ments. Typically, the Fermi energy is displaced by about 
300 meV in the experiment of Ref. jl^]. Finally, to de- 
scribe semi-conducting or insulating SWNTs with a band 
gap AE, a bulk mass term should be added |39|, 



H' = 



AE 



—paa 



(2.7) 
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We will put AE = in the following, as is appropriate 
for armchair nanotubes. 



B. Coulomb interaction 

Let us now examine Coulomb interactions mediated 
by an arbitrary potential U(r — r'). The interaction is 
described by the Hamiltonian 



<T<j' 

x U{r - r')^ a ,{r')^ a {r) . 



(2.8) 



Herein we have assumed that the interaction is not sensi- 
tive to the electron spin a — ±. Bound electrons can be 
incorporated in terms of a background dielectric constant 
re, but free charge carriers in nearby gates could lead to 
an effectively short-ranged potential. 

For the experiments of Ref. pl| , one has an externally 
unscreened Coulomb interaction, 



U(r 



e 2 /re 



-A 2 



4R 2 sm 2 [(y-y')/2R] + a 2 



(2.9) 



where a z ~ 3as ~ a [with the Bohr radius as = 
h 2 /me 2 — 0.529A] denotes the average distance between 
a 2p z electron and the nucleus, i.e., the "thicknes s" o f the 
graphite sheet. The dielectric constant in Eq. ( |2.9|) can 
be estimated from the following elementary considera- 
tion. In a strictly ID system of length L with interaction 
potential u(x — x') = e 2 / k\x — x'\, the charging energy is 

1 f L 

E r = — - / dxdx' u(x — x') 
2L 2 J V ' 

~ e 2 \n{L/R)/nL . 

The experimental value E c = 2.6 mcV [|l8| for L = 3/im 
leads to re sa 1.4. The theoretical estimate re w 2.4 for 
bulk graphite EJ is of the same order of magnitude. 
Our estimate also includes the effect of the insulating 
substrate. 

can now be reduced to a ID in- 



The interaction (2 



teraction by inserting the expansion (2.5) for the electron 
field operator. The result is 



pp'aa' Q'lQ2Q3Q4 

fct 



dxdx' 



V(l } (x-x') (2.10) 



with the ID interaction potentials 

V^ } (x-x') = J*** dydy' <p; ai (x,y)<p;, a2 (x',y') (2. 
x U(x - x', y - y' + pdS p -p>) <p p , a3 (x' , y')<p pa4 (x, y) . 



These potentials only depend on x — x' and on the ID 
fermion quantum numbers. For interactions involving 
different sublattices for r and r' in Eq. (2.8), i.e., p ^ p' , 
one needs to take into account the shift vector d = (0, d), 
see Eq. (Eml). 



To simplify the resulting ID interaction ( 2.10| ), we now 
exploit momentum conservation. Since an additional 
gate voltage tunes the average electron density on the 
SWNT, we assume that no lattice instabilities due to 
Umklapp scattering are present, i.e., one stays off half- 
filling. This assumption does not limit the applicability 
of our approach, but simplifies matters considerably. For 
a recent study of Umklapp processes in SWNTs, see Ref. 
|l2"f . We then employ the Fermi point quantum numbers 
a>i to classify the allowed processes, see Fig. ||. First, 
there are "forward scattering" (oFS) processes, where 
cti = «4 and «2 = ct3- Second, we have "backscattering" 
(ctBS), with a.\ = — oi2 = c*3 — —04. Notice that the 
above classification differs from the conventional one for 
the standard two-chain problem . The latter is based 
on the right- and left-moving indices (r = ± or p = ±) 
rather than on the different Fermi points of the nanotube 
(a = ±). 

C. Forward scattering 

Let us start with aFS, where a\ = «4 and ct2 — (X3. 
The corresponding diagram is shown in Fig. 0. We first 
define the potential 



This gives from Eqs. (|2.1l[ ) and (|J) the {aj- 
independent forward scattering interaction potential 



SV p (x) 



(2.13) 



with the correction term 



r27rR dvdv' 

SV P (x)= j 7 JLl-lU(x,y - y' +pd) ~ U(x,y - y>)} , 



{2irR) 2 



(2.14) 



which is only present if r and r' are located on different 
sublattices. Thereby important information about the 
discrete nature of the graphite network is retained de- 
spite the low-energy continuum approximation employed 
in our formulation. The correction SV P is a direct mea- 
sure of the difference between intra- and inter-sublattice 
interactions. Because o f the periodicity of the y-integrals, 
Taylor expanding Eq. ( 2.14 ) in powers of d implies that 
the correction 5V p (x) — 0. Since Vo(x) treats both sub- 
^) lattices on equal footing, the resulting aFS interaction 
couples only the total electron densities, 
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H aFS = \ J dxdx' p(x)V (x 



where the ID density is defined as 



x')p{x') , (2.15) 



paa 



(2.16) 



This density does not contain "fast" terms arising from a 
mixture of the different (paa)-type fermions. Such terms 
turn out to be cr ucial for an understanding of the Fricdcl 
oscillation in Sec. 
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For t he un screened Coulomb interaction (2.6), the po- 
tential ( 2.12 ) becomes 



V (x) = 



2e 2 



2R 



nn^/al + x 2 + 4R 2 I ^Ja 2 + x 2 + 4R 2 



(2.17) 



with the complete elliptic integral K(z) of the first kind 
IH. For x > R, one finds V (x) ~ 1/x again. The 
Fourier transform Vo(k) = Vq(— k) can be expressed in 
terms of the modified Bessel function Kq(z) [EgJ , 



Vo(k) 



4e^ 

K7T 



tt/2 



dip Kq k 



4R 2 sm 2 (p) . (2.18) 



Using asymptotic properties of K o (z), we then find the 
long- wavelength form valid for kR <^ 1, 



2fi 2 



\lnkR\ + c ) , 



(2.19) 



where c = —7 + (7r/2) In 2 ~ 0.51 with Euler's constant 
7. The logarithmic singularity at small wavevectors is a 
direct consequence of the long-range tail of V (x) . 

If one starts instead from an on-site Hubbard interac- 
tion glU, 



U(r - r') = US a ,- a ,S pp ,S(r - r') , 



(2.20) 



one finds from Eq. (2.12) the ID forward scattering po- 
tential 



V (x) 



U 
2^R 



5 (T - (T i8(x - x') 



(2.21) 



The characteristic 1/R scaling of the effective ID inter- 
action has also been found in Ref. [42| . 

For \x — x'\ 3> a, our continuum calculation leading to 
8V p {x — x') ~ is justified. However, for \x — x'\ < a, 
one has to be more careful. Here an additional aFS 



term beyond Eq. (2.15) arises due to the hard core of 
the Coulomb interaction. To study this term, we put 
x = x' and evaluate SV p (0) directly on the wrapped 
graphite network shown in Fig. [| For simplicity, we 
again consider an armchair SWNT, albeit the same re- 
sult [see Eq. ( 2.23| ) below] is found for other nanotube 
geometries. We now write 



5V p (x 



-2f p S(x 



(2.22) 



and since f p — f— p — f [see below], we obtain the addi- 
tional aFS contribution in the form 



H 



(i) 

aFS 



-/ 



pact' ' aa' 



(2.23) 



The coupling constant / > depends on the nanotube 
geometry only and is evaluated for the (n, n) case next. 

To calculate SV P (0) on the wrapped graphite lattice, we 
start from the microscopic arrangement of carbon atoms 
around the waist of the armchair SWNT, see Fig. |. On 
the length scale \y — y'\ < 2nR, it is always justified to use 
the externally unscreened Coulomb interaction potential 
fl2.9p . Denoting the locations of carbon atoms on sublat- 
ticc p = + (the circles in Fig. ||) as yk — 2n(k/n)R = 3kd 
with k = 1, . . . ,n, and discretizing Eq. (2.14), we obtain 
from Eq. ( |E2§ ) 

n 

2/ p /a = (3d/27ri?) 2 M°>W - Vk) 
-U(0,yi-y k +pd)} . 

The summation extends over the p = + sites only, and 
inserting Eq. (|2.9[) yields with R = y/3na/2ir, 



(2.24) 



Ik \fsin 2 [{l - k)n/n] + {a z /2Rf 

1 I 

^/sin 2 [(; - k +p/3)Tr/n] + (a z /2R) 2 

The singular contributions are picked up from I = k and 
I = n — k and yield the p-independent contribution 



f/a = 7/ e 2 /R , 
with the order-of-magnitude estimate 



(2.25) 



7/ 



\/3a 
2nna z 



y/1 + a 2 /3a 2 z 



0.05 . 



It can be ch ecked easily that the remaining terms in the 
summation ( 2.24 ) are neg ligibl e against the singular con- 
tribution leading to Eq. ( 2.25 ). The prefactor 7/ is in- 
dependent of n such that / ~ 1/n. Parenthetically, in 
the language of the Hubbard-like models employed in 



Refs. 



], we have / /a = U — V , where U 



is the on-site and V the nearest-neighbor Coulomb inter- 
action. According to Eq. fl2.25|) , this difference is small 
compared to U. Using a strictly on-site Hubbard model 
(V = 0) is therefore never realistic, even in the presence 
of close-by screening gates. 
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D. Backward scattering 



H = H +H 



(o) 

qFS 



H 



(i) 



aFS 



H. 



aBS 



(3.1) 



Let us now discuss the aBS contributions depicted 
in Fig. ||. For backscatterin g pro cesses, we have Oi = 
— ai = Q!3 = — 04 = a in Eq. (2.10), and from Eq. (2.11) 
the corresponding ID interaction potential for aBS pro- 
cesses reads 



VP p '(x - x 1 ) = e 2to M*-*') V pp' s ( x _ x >) 



(2.26) 



Because of the rapidly oscillating phase factor, the only 
non- vanishing contribution to H a Bg comes from \x— x'\ < 
a, i.e., we can effectively take a local interaction. Accord- 
ing to Eq. ( 2.10D , the aBS contribution must then be of 
the general form 



H, 



aBS 



dx (P+PP'b'^Ua^t'-aa'^p'aa'^p-ac 
pp' aerer 



The term corresponding to b' is irrelevant fl47| and is 
omitted in what follows. The effective coupling constant 
b > is then determined by the Fourier transform of 
Vq{x) at k = 2kp, 



(2.27) 



b = V {2k F ) , 
and the backscattering contribution is 

H aB S = ^ J dx ^Ua^l-aa'^p'aa'^p-aa ■ 

pp'aao-' 

(2.28) 

For the unscreened interaction ( |2.9| ), using 2k fR = 
4n/-\/3 3> 1 and asymptotic properties of the Bessel func- 



tion, we obtain from Eq. (2.1 



b/a = -y b e 2 /R , 
with the order-of-magnitude estimate 



(2.29) 



lb 



2tt 2 k 



0.1 



The prefactor 7b is independent of n, and thus b ~ 1/n. 
This calculation for unscreened interactions predicts b « 
/. If the Coulomb interaction is externally screened, how- 
ever, a qualitatively different situation can arise. Now / 
is still given by Eq. (2.25), but b can become significantly 
larger, see Eq. (2.27). Therefore it is possible to have 
b / in the presence of screening gates. 



III. BOSONIZATION 

According to Sec. [n| the low-energy theory of arm- 
chair SWNTs away from lattice commensurabilities is 
described by the Hamiltonian 



This model is equivalent to two spin-i fcrmion chains 
coupled in a rather special way by the interactions, but 
without interchain single-particle hopping. As described 
below, it is in that respect that our theory differs from the 
standard two-chain problem of coupled Luttingcr liquids 
investigated, e.g., in Refs. |4§|-|5~|l . The standard two- 
chain model is usually derived by coupling two Hubbard 
chains via a transverse hopping matrix element |5q| , and 
has been used to study, e.g., the stability of Luttinger liq- 
uid behavior with respect to interchain coupling [p r 6|-^8| . 
In contrast, the modified model (3.1) describes the prop- 
erties of a SWNT. Our solution of this model proceeds 
in fact quite similar to the standard two-chain case. In 
particular, we now follow the bosonization route taken 
in Refs. p9|-pl| and especially in Ref. |5^] . For a general 
review of bosonization, see Ref. J2j| . 

In order to proceed, we need to bring the non- 



interacting Hamiltonian (2.6) into a standard form of the 
ID Dirac model. This is accomplished by switching to 
right- and left-movers (r = ±) which are linear combi- 
nations of the sublattice states p = ±. In particular, we 
have 



(3.2) 



where the unitary operator U fulfills U^a y U — a z . Since 
the indices of U are p — ± and r — ±, the 2x2 matrix 
representation U = (l/\/2)[(l, «), (1,— i)] is easily found. 
The non-interacting Hamiltonian (2.6) then reads 



H = -iv J dx i\)\ a(J d x ^) r 



(3.3) 



Furthermore, the ID density operator p paa -(x) entering 
Eq. (2.16|) becomes 



r=± 



(3.4) 



where p ra <j( x) — i^l n/T (x)ip rnn (x). The ID density de- 
fined in Eq. ( 2.16|) is therefore equivalently expressed as 



(3.5) 



In this representation, it is convenient to apply the 
bosonization formula |E3l, 



exp {iqprx + ikpax + iip ra cr} ■ (3-6) 



For simplicity, we have incorporated the spatial depen- 
dence due to the Bloch functions (2.2) into the ID 
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fermion operator. The chiral bosonic phase fields obey 
the commutator algebra 

[ip ra a(x),(Pr'a'a'(x')]_ = -iltr8 rrl 8 aa , 8 aa > Sgn(x - x') . 

(3.7) 

From these relations, the fermion operator rp raa (x) has 
indeed the correct anticommutator algebra on a given 
branch (rata). To enforce anticommutator relations be- 
tween different branches, we need the Majorana fermions 
(Klein factors) i] raa . — r)l aa .. They obey 



\f}ra.a , Tjr'a'a' } _|_ 2$ rr > S aa ' 8aa' 



(3.8) 



We note that the usual zero modes p3| have to be incor- 
porated in the fields tp raa as they are not contained in 
the Majorana fermions. However, they matter only i f on e 
is interested in finite-size effects. Finally, from Eq. (3.6), 
the density p ra a now takes the form 

If 



Praa (^) 



2?r 2tt 



dx^Praa (^) 



The wave vector qf, which must be carefully distin- 
guished from the Fermi vector kp, is related to the band 
filling. By varying a suitable gate voltage plf , one can 
easily adjust the Fermi energy and hence the band fill- 
ing (see Fig. |l|). Relative to the unbiased half- filled case 
Ep = 0, an average excess density Sp gives rise to a 
non-zero qp = irSp/A and hence Ep — vqF. We con- 
sider \qp | <C kp in the following, since otherwise the 
low-energy continuum approximation underlying our ap- 
proach might break down. 

At this stage, it is natural to introduce the standard 
linear combinations 9j${x) and their dual fields <pjs(x) 
subject to the algebra 

[<j) :j s(x),9 jl s'(x')} = -(i/2)5jf5 S s' sgn(x - x') . (3.9) 

The bosonic phase fields 9j$ (x) for the total and relative 
(5 = ±) charge (j = c) and spin (j = s) channels are 
constructed as 



1 



7 c+ 



40? 
1 

1 



^ ra ip ra(7 , 

rota 



7 s+ 



raa 
1 X - 



40? 



raa 



ra lfi r 



raa tp racr 



Their dual fields (j)js are defined similarly by omitting 
the r factor in these relations. The back-transformation 
reads 

fraa = "^r{^+ + r #e+ + 01 <fi c - + ra6 c - (3.10) 



+ a <fi s+ + ra 6 S+ + aa < 



} 



such that Eq. (3.7) is recovered from Eq. (3.9). The total 
electron density relative to the half-filled situation reads 
from Eq. (§13) 



p(x) = Aq F /-K H =d x 9 c+ (x) , 

V7T 



(3.11) 



and the continuity equation then yields the current, 
2e 



i = ^d t e c+ (x,t) 



(3.12) 



which can be evaluated, say, at x = 0. 

In order to arrive at the bosonized form of the model 
( |3.lD , we need to specify the Majorana fermions r\ raa in 
Eq. ( |3.6[ ). Since the Hamiltonian contains only the spin- 
conserving products 



A 



(3.13) 



these can be represented using standard Pauli matrices 
p^ , p9[ . Besides A ++ = 1, we choose 

= iaa x , A |_ = iraa z , A = —ira y . (3.14) 



To show that Eq. ( 3.14 ) leads to a valid representation for 
the Klein factors, we have to check all possible products 
of A±± with each other. For instance, the product 



— A — \- — 7]racr11r—a(T'r]ra(TT]—7 



(3.15) 



gives (iaa x )(iraa z ) = ira y according to Eq. ( 3.14 ). On 
the other hand, using the anticommutator relation (3.S) 
gives for Eq. (3.15 ) —r\r~aoT}-ra<j, which is again ira y 
from Eq. (3.14). Similarly, all other products are con- 
sistent with the algebra (3.8), and one can indeed use 
the representation (3.13) with Eq. (3.14). Of course, our 
choice (3.14) is not unique, and one can find other possi- 
bilities. However, final expressions for physically observ- 
able quantities do not depend on this choice. 

The bosonized expressions for the various terms in 
Eq. |H]) read @ 



dx 



K jS {d x <t> o5 f +Kj s \d x 6 jS ) 



H 



is 

i% = - / dxdx ' d x 6 c+ (x)V (x - x')d x ,6 c+ {x') 



f 



■ J dx [- cos(\/47r 6> c -) cos(V47r 6> 5 ._) 



(i) 

,ll/s ~ (™) 2 

cos(VAn 9 C -) cos(V47r 9 S+ ) 



H 



+ cos(v47r6> s+ )cos(V47r a _)] 

H a BS = t — s~ / dx [cos(a/47t0 c -) cos(v / 47r6' s _) 
{way J 

+ cos(\/47r6' c _) cos(V47T0 ; ,._) 
+ cos(\/47r# s _) cos(v / 47r</> s -)] ■ 



(3.16) 
(3.17) 

(3.18) 
(3.19) 
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The correct signs for the various terms in Eqs. (3.18) and 
( |3.19 ) are crucial and necessitate a correct choice for the 
Majorana fermion products, see Eq. (3.14). Although 
bosonization of Eq. (3.3) gives Kjs = 1 in Eq. ( 3.16 ), 
interactions renormalize these parameters. In particular, 
in the long- wavelength limit, H^ s can be incorporated 
into Hq by putting 



K c+ =K= h + 4y (k~Q)/nvX V2 < 1 , (3.20) 

while for all other channels [(jS) ^ (c+)], the coupling 
constant / gives rise to the renormalization 



Kj S = 1 + f/TTV > 1 



(3.21) 



The corresponding renormalization of K c+ due to / ^ 
can be neglected against the large effect of E q. (3.17 ). 
The velocities of the different modes in Eq. (3.16) are 
then given by 



Vjs = v/Kjs . 



(3.22) 



Clearly, the charged (c+) mode propagates with signif- 
icantly higher velocity than the three neutral modes. 
There is a further renormalization of the values given 
by Eq. ( [3.22 ) due to the aBS contribution. However, 
this effect is very small and does not affect the power 
laws discussed below. We shal l the refore neglect velocity 
renormalizations beyond Eq. ( 3.22) in t he following. 

For the long-ran ged in teraction ( [2.19 ), the logarithmic 
singularity in Eq. ( |3.20 ) has a natural infrared cutoff at 
k = 2it/L due to the finite length of the nanotube. For 
temperatures ksT ^> hv/L, thermal effects provide a 
higher cut-off. Considering low temperatures, ksT <C 
hv/L, we then obtain 



K = 1 1 



8e 2 l - 1/2 

— ^— [ln(L/27ri?) + 0.51] 

7TKHV 



(3.23) 



which predicts K — > for L — > oo. Since Tic/e 2 ~ 137, 
we get with v = 8 x 10 5 m/sec the estimate e 2 /Tiv = 
(e 2 /hc)(c/v) w 2.7, and therefore K ~ 0.18 for the L = 
3fim tube of Ref. Jig ], Quite generally, the parameter K 
give s the appropriate measure of the correlation strength 
j22"| , ^3f . The rather small value of K found here implies 
that a single-wall armchair nanotube is a very strongly 
correlated system, which should display pronounced non- 
Fermi liquid effects. 



IV. RENORMALIZATION GROUP ANALYSIS 



The Hamiltonian given by Eqs. (fj(|-(fl§) does not 
allow for an exact solution. To proceed, we investigate 
the nonlinear terms associated with the coupling con- 
stants / and b by using the perturbative RG method 



[pl|| . Since there is no renormalization in the charged 
(c+) sector, in the following (jS) includes only the com- 
binations (c— ), (s+), and (s— ). There are eight scaling 
operators perturbing the critical b = f = model that 
need to be taken into account. They are given by 



Vi = 


1 

— cos(v 

7T 




) cos( i 




(4.1) 


v 2 = 


1 

— cos(v 

7T 




) cos( i 




(4.2) 


v 3 = 


1 

— cos(v 

7T 




) cos(- 




(4.3) 




1 

— cos(v 

7T 




) cos( i 




(4.4) 


v 5 = 


1 

— cos(v 

7T 




) cos( i 




(4.5) 


v jS = 




+ (d x t 


^) 2 ] > 


(4.6) 



with associated scaling fields (coupling constants) gi(£) 
(where i = 1, . . . , 5) and gjs{t). Here dl = — din ui c de- 
notes the standard RG flow parameter, i.e., one decreases 
the high-energy bandwidth cutoff uj c and compensates 
this decrease by adjusting the coupling constants. Within 
the reach of the perturbative RG, the coupling constants 
gjS are related to the Kj$ parameters by Kjs = 1 — 2gjs- 
The initial values of the coupling constants are 



Si(0) = (&-/)/™, 
52 (0) = - 53 (0) = 2g jS (0) 
5 4 (0) = b/irv , 
.95(0) - . 



-f/iru , 



The operator V5 is not present in the original Hamilto- 
nian but will be generated during the renormalization 
process. Furthermore, we have omitted the operator 

V' = -cos(\/4^6> s -)cos(\/47r0 s _) , 



which is present in Eq. ( 3.19 ). This operator stays ex- 
actly marginal in all orders and thus decouples com- 
pletely from the problem. In fact, by using the Majorana 
refermionization procedure in Sec. [y], this term is seen to 
vanish. 

It is straightforward to derive the second-order RG 
equations, e.g., by using poor man's scaling or, more el- 
egantly, the operator product expansion technique J6l] . 
The resulting equations read 

dgi 
dl 

dm 
dl 

dg 3 
dl 
dgi 
dl 



5253 + gi (g c - 


+ 9s-) , 


(4.7) 


5153 - 5455 + 


52 (5c- +5s+) . 


(4.8) 


5152 + g3(g s + 


+ 5,-) , 


(4.9) 


5255 + 54 (ffc- 


- 9s-) , 


(4.10) 
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dg 5 . , 

—j£ = ~929i + 95{g s + - 
dg c - 1/2 , 2 i 2\ 

^ r = 2(.g 1 +.9 2 + 3 4 ) 



9s-) , 



dg. 



s + 



dl 

dg s - 
dl 



\{gl + gl + gl) , 
\{g\ + gl-gl-gi) 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



We shall first study these equations in two limiting cases, 
namely either b — or / = 0. 

Let us start with the case b — 0, where the coupling 
constants 34 and 35 stay identically zero. Putting 

G = -.91 = -.92 = .93 

with the initial condition G(0) = f/nv > 0, and 

G = 1g c - = 2.9s- = 2g s+ 

with G(0) = —f/irv < 0, the RG equations read 



dG 

~dJ 



= -G 2 



GG , 



dG 

~dJ 



= 2G 1 



This permits the solution G — — G obeying dG/dl = 
— 2G 2 . Since G(0) > 0, all coupling constants flow to 
zero as I — > 00, 



G(£) = 



G(0) 



1 + 2G(0)£ 



Therefore, for 6 = 0, the contribution -ff^pg in Eq. (3. IS) 
is marginally irrelevant. 

Next we turn to the case / = 0, where one can imme- 
diately put 32 = 33 = 35 = ,9s+ = 0. Furthermore, we 
are free to choose the isotropic solution 

.9 = ,9i = 34 , 



with the initial condition g(Q) = b/nv > 0. Then one can 
also put g s - = and is left with the Kosterlitz-Thouless 
equations 



dg 



gg c - 



dg c - 
d£ 



(4.15) 



Therefore both g and g c - will flow to strong coupling, 
and the backscattering part i? Q BS is marginally relevant. 
The physical picture emerging for the special case / = 
is elucidated in Sec. From this analysis and the nu- 



merical solution of the general RG equations (4.7)— (4.14) 
for / 7^ 0, it is apparent that b is always a marginally rel- 
evant coupling constant. 

The next question to be addressed is whether the ini- 
tially irrelevant H^fs becomes relevant near the emerg- 
ing strong coupling fixed point. To study this issue, let us 
assume b 3> / and linearize the RG equations in the ini- 
tially small coupling constants 32,^3,^5, and gj$. Then 



9 S + decouples, and we introduce the new coupling con- 
stants 

.9 = (31 + 54)72 , g = gi - gi , 



with initial values g(0) — bjixv and g(0) = —f/irv. 
Thereby we recover the Kosterlitz-Thouless equations 
(4.15), so that the above strong-coupling behavior of b is 
basically unaltered, i.e., the coupling constants g and g c _ 
flow to strong coupling, g — > g* > and K c ^ — > K*_ < 1. 
The remaining RG equations read 



2gg s 



dg 
dl 
dg2 
d£ 

dg3 _ dg 5 
d£ " d£ 
dg s 



3(33+35) , 



"332 



d£ 



gg 



Therefore, g(£) also reaches the strong coupling regime, 
g — ► g* < 0, and similarly, K s - — * K*_ > 1. Fur- 
thermore, it follows readily that 32,33 and 35 flow to 
strong coupling (with 33 — 35 staying constant). We con- 
clude that terms associated with the coupling constant 
/ will become relevant near the strong-coupling point 
b — > b* , which therefore represents only an intermediate 
fixed point. 

This conclusion remains valid if the assumption b 3> / 
is relaxed. In Fi g. fi| , t he nu merical integration of the RG 
flow equations (4.7)-( 1.14 ) is shown for realistic values 
of the coupling constants b and /. Clearly, all coupling 
constants go to the strong-coupling regime, albeit 34 asso- 
ciated with the backscattering contribution is dominant. 
Only at sufficiently low energy scales, i.e., for large values 
of the flow parameter £, the nonlinear forward scattering 
term -ff^pg becomes important. Finally, in Fig. [|we show 
the RG flow for different initial conditions corresponding 
to / > b. The coupling constants 32 and 33 associated 
with the forward scattering contribution are initially ir- 
relevant. However, they eventually become relevant and 
flow to strong coupling as £ — * 00. In the end, all coupling 
constants flow into the strong-coupling regime again. 



V. MAJORANA REFERMIONIZATION 

We start this section by presenting a solution for the 
special case / = 0. This solution proceeds by Majorana 
refermionization and is similar to the one discussed by 
Schulz in the context of the standard two-chain problem 
j52j. For / = 0, the nonlinearity affects only the relative 
(5 = — ) channels, for which the bosonizcd Hamiltonian 
reads from Eqs. ( ft. If ) and ( [3.19 ) 
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H(c-,s-) = lJdxJ2 [(d^j-f + (9Jj-) 2 ] (5.1) 

j = C,S 

+ - — — / dx [cos(V47r(9 c _) cos(V47r# s _) 
{ira) z J 



4- cos(v 47r 9 C J) cos(v 47r </> s _) 
+ cos(V47r0 s -)cos(v / 47r<^ s _)] . 

Let us define new effective fermion operators for the rela- 
tive charge and spin (J — c, s) channels. Their right- and 
left-moving components (p = ± = R, L) can be written 
in terms of the bosonic phase fields, 



Then we have 
cos[\/47r0j-] = 



jhp_ 

J2na 



exp{—i^(p6j- + <t>j-)} 



= -T] jR r) jL ira (^) R ^] L - ^jL^jr 



(5.2) 



Next we express these two (complex) Dirac fermions in 
terms of four (real) Majorana fermions £,j P (x). Here j = 
1,2 corresponds to the (s— ) channel, and j — 3,4 to 
(c-), 



The Majorana fermion operators obey the algebra 
[€jp{x),£j> P '(x')]+ = 5jf5 PV ,8(x - x) , 



(5.3) 



and by using Eq. (5.2) we find 

cos[\/47r6> s _] = -ina (£ir£il + &r&l) 
cos[V47T^ s -] = -«7ra {CirS.il - 
coslv^^c-] = JTra (£, 3 r£.3L + UrUl) , 



(5.4) 



where the Klein factors have been chosen as rj cR rj cL = 
—rj sR r) sL = —i. This choice is dictated by the condi- 
tion that the interaction term is marginally relevant, see 
below. 

As already mentioned, one can see from this represen- 
tation that the operator 



COs[v47T0 s _] C0S[V 47T0 S _] 

can effectively be put to zero. To be precise, upon point 
splitting [p3|]6"T|| , it yields highly irrelevant operators in- 
volving spatial derivatives o f the Majorana fields. The 
corresponding term in Eq. (5.1) can thus be omitted. 
Refermionization then yields [we put v = 1 in the in- 
termediate steps] 



3=1 



H (c-, O = ~\ X! / dx : s//.'' v 's//.' - <><!;, i. ) 



+ 2bjdx i&R&L + UrUl) tiR&L • (5.5) 

In order to better understand the above model, it is 
instructive to rewrite it in terms of the following current 
operators, 



I R (x) = i£x R (x)£ 3R (x) , 
7 fl( x ) = -i£iR(x)Zm(x) , 

I R (x) = i£,3R(x)^4R(x) , 



(5.6) 



and analogously for the left-movers. These currents obey 
the SU(2) level-2 Kac-Moody algebra ||, 

[I a (x), I\y)]- = ie abc 5{x - y)I c (x) + ^<5'(x - y) . 

Z7T 

(5.7) 

In terms of the operators ( |5.6| ), the interacting part of 
the Hamiltonian (|5.5|) takes the form 



dx {g ± [r R (x)IZ(x) + I y R {x)Ii{x)] + g z I R (x)It(x)} 



The transverse coupling is g± = 2b, while the longitu- 
dinal coupling g z is initially absent but will be gener- 
ated under the renormalization. We thus can identify 



the model (5.5) as the SU2{2) anisotropic Wess-Zumino- 
Witten model, see, e.g., Ref. |2j|. From the struc- 
ture of the commutation relations, it immediately follows 
that the coupling constants (g±,g z ) obey the Kosterlitz- 
Thouless equations 



dg± _ _L_ 

di 2n 

Using the operator identity 
1 



dg z 
di 



cos[v47r#(a; 



1 2 

Z7T 



■ cos[vT&7T0(a;)] 



(5.8) 



-(d x 9f + const. 

7T 



(5.9) 



we conclude that the longitudinal coupling corresponds 
to a renormalization of K c ^, with g z = 2ng c ^. Using 
this identification we recover the flow equations (4.15) 
and highlight the underlying reason for their Kosterlitz- 
Thouless structure. 

Of course, for b > 0, the model flows to strong cou- 
pling, rendering the Majorana fields £i, £3 and £4 mas- 
sive. Notice, however, that the Majorana fermion £ 2 
stays massless, so that the (s— ) sector carries one massive 
(j = 1) and one massless branch. This behavior is due to 
the symmetric appearance of the dual fields # s _ and </> s _ 



(self-duality) in Eq. (5.1). By virtue of the Heisenberg 
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uncertainty relation, it is impossible to completely pin a 
self-dual field. However, the (c— ) sector is fully massive. 

The masses of the massive branches can now safely be 
calculated by Majorana mean-field theory. We mention 
in passing that a standard self-consistent harmonic treat- 
ment p| does not apply in this case since we have to deal 
with a marginally relevant perturbation. In the mean- 



field approximation, the interaction term in Eq. (5.5) is 
written as 



2ib J dx (ci [£ 3 fl6z, +Ur£±l] + 2c 3 CliiCii) 



(5.10) 



with the ground-state averages 



Cl 
C 3 



j(£ifl£iL> > 



4L 



At this point it is straightforward to diagonalize the 
mean-field Hamiltonian. This yields the self-consistency 
relations 



ci 



26c, 



bci 



ln(w c /46c 3 ) 



c 3 = — ln(w c /26ci) , 

where ui c = 7.4 eV is the bandwidth of the ir electrons 
j|. The solution to these relations is c 3 ~ ci/a/2 with 

ci — 2^ exp[-7r/V2&] . 



Comparing Eq. (5.10) with the usual mass term 
— im J dx^n^L yields the masses nij of the massive Ma- 
jorana fermion branches j = 1,3,4. Apart from a factor 
\/2, these masses are all equal and given by 



nib = exp[— t:v /V2b] 



(5.11) 



Since the coupling constant b appears in the exponent 
of an exponentially small quantity, it is quite difficult to 
estimate the value of mj. Taking the value (2.2E) for b 
yields an order-of-magnitude estimate for the associated 
temperature scale, Tf, « 0.1 mK. For well-screened inter- 
actions, however, Tb can be orders of magnitude larger. 

To proceed, since the (c— ) sector is fully massive, we 
now effectively put Q c - = in Eq. (5.1). Other degener- 
ate solutions give the same physical results. Then we are 
confronted with the interaction Hamiltonian 



2na 



J dx ^cos[V47T# s _] + cos[\/47r0 s _]^ 



Self-duality implies that sin[y / 7r# s _] and sin[Y / 7r(/> s _] must 
have the same scaling dimension. Similarly, the scaling 
dimensions of cos[-\/7r# s _] and cos[^/tt(/) s ^] must coincide. 
Although the fields # s _ and <fi s - are not pinned, they still 
have a tendency to be pinned. Since nib > 0, this implies 
that the values of sin[y/ir9 s -] and sin [y^^s-] are close to 



1, while the corresponding cos-operators approach zero. 
Therefore the sin-operators will be correlated stronger 
than the cos-operators. 

These simple arguments can be made quantitative as 
follows. The (s— ) sector is described by two decoupled 
Hamiltonians, H = //i[£i] + #2 [£2], for the Majorana 
fermions £1 and £2, 

HjM = ~\ j ^ i&RdvtjR - fcifla&z) (5.12) 

- inij I dx i, /,■£,/. , 



where the masses are given by mi = nib and nii = 0. 
We now exploit the well-known correspondence between 
the 2D classical Ising model and ID Majorana fermions 
p3| , |62[ . The fermion mass nij sets the relevant energy 
scale, where the correlation length of the related Ising 
model should be proportional to nij . Therefore the Ising 
model related to £2 will be critical (T — T c ), while the 
Ising model for £1 is above criticality (since mi > 0, we 
have T > T c ). Denoting the standard order operators for 
the two Ising models as <j\ 2 , and disorder operators as 
fj.1,2, the fusion rules [{23|]62t yield the correspondence 



cos[y/n9 s - 

cos[\/ir4>s- 
sinfv^s- 
sin[v7n/> s _ 



= 01M2 
= C1C2 
= Atl02 

= 



(5.13) 



These fusion rules state that for two given Ising models, 
products of the order/disorder operators determine the 
Majorana field operators, which in turn are composed out 
of the exponentials cxp[±iy/Tr9 s -] and exp[±«y / 7r0 s _]. 
The relation ( 5.1 3| ) is a valid representation of the fusion 
rules, see Ref. ]23j|. 

Since the first Ising model is above criticality, the aver- 
age of the order operator is zero, (<7i) = 0. However, the 
disorder operator /xi has a finite average value, and cor- 
relation functions of cos[y / 7r^ s _] and cos[\/7r</) s _] decay 
exponentially. From the exact solution of the 2D Ising 
model 62 one then immediately obtains 



(sin[v / 7T^ s -(a;)] sin[\/7r0 s _(x')]) ~ \x — x 



'1-1/4 



(5.14) 



with the same result for the sin[Y / 7T0 s _] correlator. The 
sin-operators therefore contain the Ising operator /ii from 
the off-critical model which has a finite average. The 
other Ising operator then comes from a critical model 
(m 2 = 0) and leads to the scaling dimension n = 1/8, 
see Eq. (5.14). In contrast, the cos-operators contain 
the Ising operator o\ with zero expectation value. Note 
that for b = all these operators have scaling dimension 
77 = 1/4. The "halved" scaling dimension n — 1/8 found 
for b > can be traced to the fact that one Majorana 
fermion becomes massive while the other stays massless. 
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In the (c— ) sector, correlations of cos[y/T:9 c -] show long- 
range order while all other operators lead to exponential 
decay. 

The results obtained so far by Majorana refermion- 
ization apply only to the intermediate fixed point char- 
acterized by / = 0. As discussed in Sec. [V, while 
-ff^pg is irrelevant around the non-interacting fixed 
point, it becomes relevant near this intermediate strong- 



coupling point. The term ~ cos(v4^^s+) cos(v47r6' s _) 



m Eq. ( |l8| ) stays marginal, but the two other terms be- 
come relevant with scaling dimension rj = 1 . This renders 
the (s+) channel massive. Regarding the (s— ) channel, 
the most important contribution due to H^p S comes from 



H 



(i) 

qFS 



/ 



(na)'- 



■ J dx cos[\/47r# c _] cos[\/47r6> s _] 



Noting that we have chosen # c _ = 0, we can again em- 
ploy Majorana refermionization. The only but important 
effect of this contribution then consists of a renormaliza- 
tion of the masses mi^, 

mi — > (mb ~ m f) j m 2 ~> m f i 

which breaks the self-duality in the (s— ) sector and thus 
drives the second Ising model off criticality as well. Ma- 
jorana mean-field theory yields 



m f — (f/b) m b 



(5.15) 



and the associated temperature scale Tf w (f/b)Tf,. 
Therefore both the Majorana fermion £2 and the (s+) 
field acquire the (generally small) mass m/ due to the 
forward scattering contribution. At the emerging T = 
strong-coupling fixed point, we have long-range order in 
the operators 

cos[v^r6' s+ ] , cos[^/~k9 c -] , sm[y/TT(f> s ^} , (5.16) 

with exponential decay in all other operators except those 
of the critical (c+) sector. That the first two operators 
become long-ranged is a direct consequence of the pin- 
ning condition # c _ = 0. That sin[y / 7r0 s _] exhibits long- 
range order can be easily understood from Eq. ( 5.13 ), 
since the second Ising model is now also above crit- 
icality. Therefore both (fi\) and ((12) arc finite, but 

(<Tl) = (<F2) = 0. 

This analysis for a general interaction potential pre- 
dicts that the exponents corresponding to the first (inter- 
mediate) strong-coupling point should be observable on 
temperature scales Tf < T < Tf,, with a crossover to a 
regime T < Tf dominated by the true T = fixed point. 
For long-ranged interactions, we have Tf m Tf,, and 
the intermediate fixed point and the associated crossover 
phenomenon are not observable. However, for screened 
interactions, Tf, can be significantly larger than Tf. Fi- 
nally, in the regime T > Tf,, which also includes the 



temperature range experimentally studied in Ref. [JL8j, 
the physical behavior emerging from our analysis is best 
characterized as that of a Luttinger liquid with an addi- 
tional flavor index. In this regime it is justified to neglect 
the nonlinearities associated with the coupling constants 
b and /. Nanotubes therefore constitute a realization of 
the Luttinger liquid at sufficiently high temperatures. 



VI. SUSCEPTIBILITIES AND CORRELATION 
FUNCTIONS 

With the strong-coupling solution of Sec. [v| at hand, 
we can now examine temperature-dependent susceptibil- 
ities and other experimentally accessible quantities. Due 
to the Mermin- Wagner theorem, there can be no ordered 
state in a ID system p3| ]. It is therefore customary to 
characterize the physical behavior by the most slowly de- 
caying correlation function, which in turn indicates in- 
cipient instabilities. We have investigated correlators of 
charge-density wave (CDW), spin-density wave (SDW), 
and superconducting (SC) type. The results reported in 
this section are summarized in Table |. 

One has to carefully distinguish the three temperature 
regimes discussed in Sec. [v|, since different instabilities 
emerge in different temperature ranges. Furthermore, 
we also have to distinguish the spatial oscillation period 
of the correlations. It is apparent from the dispersion 
relation in Fig. [j] that the wavelengths 



\ = ir/kp, ir/qF, 7t/(A;f ± C[f) 



(6.1) 



could occur in a doped armchair SWNT (we put qF > 
for simplicity). Without doping, only the standard 7r/fcp 
wavelength is found, which corresponds to an order op- 
erator effectively involving two ID fermions at different 
Fermi points (a = —a'). In the doped case, the much 
longer wavelength n/qF arises if the ID fermions are at 
the same Fermi point but move in opposite directions 
(r = —r'). Finally, the wavelength ir/(kp ± <?f) corre- 
sponds to the mixed situation (a — —a',r = —r'). The 
possible simultaneous occurrence of different wavelengths 
is a remarkable feature of doped nanotubes which has 
no analogue in the standard two-chain problem. This 
phenomenon has its ultimate origin in the unique band 
structure of graphite. 



A. CDW correlations 

In calculating the CDW correlations, we first need to 
find the bosonized representation of the y-integrated [this 
only amounts to a factor 2nR] density operator q(x). It 
can be derived from Eq. (|2.5|), 



q(x) 



J a 
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= J dyip* pa {x)tp p , a ,(x) ipl a<T (x)ip p , a , a (x) . (6.2) 



pp' aa'a 



There is first a "slow" component p(x) du e to p — p',a — 
a' , whose bosonized form is given in Eq. (3.11). Further- 
more, there is an intra-s ubla ttice (CDWO) order param- 
eter contributing in Eq. ( |6.2| ), 



Ocdwo{x) 



r(x) 



Since the Bloch functions for different sublattices are or- 
thogonal, one might now conclude that p 7^ p' cannot 
give a contribution to q(x). However, the corresponding 
finite (but generally small) matrix element allowing for 
a contribution due to the inter-sublattice (CDW71-) order 
parameter 

6 C DW-k{x) ~ J2tpl a o(x)^-p,±a,a(x) 
pacr 

can be generated by the interactions. Such a mechanism 
is well-known from the study of 4fcp components in the 
density operator for correlated ID fermions p3] . 

To make use of the bosoni zed version, we then employ 
the unitary transformation (3^2) and find 



Ocdwo 



~ sin[v / 7T'/ > c 
x cos[^6 s 



+ 2kFx] cos[y/n0 c -] 
I cos[\/7r</> s _] — (cos <• 



(6.3) 



sin) 



The correlation function of Ocdwo thus has the wave- 
length A = n/kp and decays exponentially for T < Th- 
in contrast, for T > TJ,, its scaling dimension is rj = 1. 
We therefore omit the CDWO mode in the sequel as it 
exhibits a fast subdominant decay for all temperatures. 

Turning to the CDW7T correlations, we find two con- 
tributions, 



61 ~ Y,(~ ir ^™<T$-r, a ,<r 

rota 

~ cos[v / 7r^c+ + 2qF%] cos[\/7r0 c _] 

x sm{^/iT0 s+ } sin[y/n9 s -] + (cos «-* sin) , 

rota 

~ cos[v / 7r#e+ + 2qFx] COS [y/n<j)c- + 2kFx] 
x cos[y/n6 s+ ] cos[y/ncj) s -} + (cos <-» sin) . 



(6.4) 



The first operator leads to a slowly oscillating correlation 
function with wavelength A = Tr/qp, while the second op- 
erator implies rapid oscillations with A = n/ (fop ± qp). 
Both operators exhibit exponential decay at the lowest 
temperatures, T < Tf. In the intermediate tempera- 
ture regime, Tf < T < O2 also leads to exponentially 



decaying correlation functions, but 0\ has the scaling di- 
mension r) = (3+2K)/8. Finally, in the high-temperature 
Luttinger liquid regime, T > the scaling dimension 
i] = (3 + K)/4 arises for both operators. For scaling 
dimension rj, the equal-time correlation function is 

{6 1 {x)6 1 {x 1 )) ~ cos[2gj(a: - a:')] \x - x'\~ 2ri , (6.6) 

and the corresponding susceptibility has the temperature 
dependence x(2<7f) ~ J 12 ^- 2 . p Qr corr elations of O2, 
instead of 2qp one has the wavevector 2{kp ± qp). Since 
only the slowly oscillating contribution due to 0\ exhibits 
power-law behavior for Tf < T < Tf,, it is favored over 
the O2 contribution by a larger prefactor for T > Tf,, at 
least for well-screened interactions. 



B. Higher-order CDW correlations 

There is also a contribution to the density correlation 
function effectively originating from squaring the above 
order parameters. Even if the original operators are irrel- 
evant, the emerging higher-order operators can become 
relevant. The analogue of this behavior in the stan- 
dard two-chain problem has been interpreted as incipient 
Wigner crystal behavior j52j , since it is characterized by 
a wavelength 7r/2fci? corresponding to the average elec- 
tronic spacing. Such an interpretation is not possible in 
our case, since the characteristic wavevector is now 4qp 
or even 8qF instead of the standard value 4/cf. 

The operator O CDW0 stays always irrelevant, but the 
squared CDW7T operators give important contributions. 
Using Eq. (5.9), we find the leading terms 



0{ ~ cos[V47r6> c+ 
+ cos[A/47r0 c _ 
60 — cos[\/47r0 c+ 



COS[V47T0 C _ 



+ Aqpx] I - cos[\/47r0 s -] 
- cos[V47r6> s+ ] j , 
+ Aqp x] jcos[%/4lr6> s+ ] 
+ -ikpx] + cos[VAtt if) s -] I 



(6.7) 



(6.8) 



Here we have kept only contributions containing the C+ 
field, since these are the only ones which could become 
relevant. Since these operators come from products of 
four ID fermion operators, the Pauli principle allows ex- 
actly three terms for 0\ and 0\, respectively. The pref- 

are 



g\ actors for the various terms in Eqs. (6.7) and (6 



non-universal and depend on the interactions. The scal- 
ing dimension of these operators is 77 = 1 + K for T > T&, 
so that they are always irrelevant in the Luttinger liquid 
regime. Furthermore, for T < T^, we observe that 0\ 
also produces only subleading contributions. In effect, 
we need to consider only 0\, which indeed becomes a 
relevant operator for T < T^. The scaling dimension in 
this temperature regime is 77 = K, leading to the Aqp 
oscillatory correlations 
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{O CDWiz {x)O cdw ^{x')) ~ cos{iq F {x - x')] \x - x'\ 



i\-2K 



(6.9) 



u CDWtt 



COS 



This turns out to be the dominant instability at T < T\, 
and K < 1/2, see Table |. 

In the same line of reasoning, there also exists an ef- 
fective contribution due to Ocdv/tt- This operator is 
relevant for extremely strong correlations, 

16^e c+ + 8q F x] , (6.10) 

with wavelength it /Aqp and scaling dimension r\ = AK for 
all temperatures. By comparing to the scaling dimension 
77 = (3 + K)/A of 61, the A = ir/Aq F CDWtt state is 
seen to represent the dominant incipient instability in the 
Luttinger liquid regime T > Tb for very strong but still 
realistic correlations, K < 1/5. At temperatures below 
T(,, the 8qp operator (3.10) leads to only subdominant 
contributions and can safely be ignored. 

C. SDW correlations 

Next we turn to spin- density wave correlations. Due 
to the underlying SU(2) spin isotropy, it is sufficient to 
study the s z correlations which can be readily evaluated. 
Similar to the CDW case, we have the order parameters 

-p,±a,a\ x ) ■ 

pacr 

Again the intra-sublattice SDWO component can be ig- 
nored as it only leads to 17 > 1. This is apparent from 
the bosonized representation 

raa 

~ cosfv^c- + 2JtFx] COs[v / 7T^c-] (6-11) 
x cos[^/tt9 s -} sin[v / 7T0 s _] — (cos <-> sin) . 

The important part comes from the SDW7T order opera- 
tor 

— p : ±a,a " 

paa 

Rotation to the right/left-moving basis and subsequent 
bosonization gives the two operators 

raa 

~ sin[ v / ^6' c+ + 2q F x] cos[V^^ c -] (6.12) 
x cos[t/tt0 s +] sm[y/n9 s ^} + (cos <-* sin) , 
O b ~ 5Z(-im)^ ao .^_ r _ a)CT 

raa 

~ cos[^/tt9 c+ + 2qpx] sin[v / 7r</' c - + 2kpx] (6.13) 
x cos[^/tt9 s +] sin[v / 7T0 s _] + (cos <-» sin) . 



Both operators lead to exponentially decaying correla- 
tions at T < Tf. Furthermore, in the regime Tt < T < 
Tfc, O a has the same scaling dime nsio n r\ = (3 + 2K)/8 
as the CDW-7T operator 0\ in Eq. (6.4), with exponential 
decay in Ob- Finally, for T > both O a ,b have again 
the same scaling dimension 77 = (3 + K)/A as 0\2- We 
then arrive at exactly the same power laws for SDW and 
CDW correlations. 

These power laws represent the most slowly decaying 
correlations of the system in the wide temperature range 
T > Tr, provided the correlations are not too strong, see 
Table |. To discriminate among CDW and SDW corre- 
lations, we have to study prefactors of the power laws. 
From Eqs. f[J) and ( |6.12| ), since cos{y/n9 s +] appears in 
O a instead of sin[y / 7r(9, s -|-] for Oi, and because at temper- 
atures T < Tf we have the pinning condition # s + = 0, 
see Eq. (5.16), the SDW amplitude should be larger than 
the CDW amplitude also for T > Tf. The magnetic 
SDW correlations then represent the dominant instabil- 
ity of the nanotubc in the temperature regime T > Tf 
for moderate correlation strength. Remarkably, they are 
characterized by the simultaneous presence of the wave- 
lengths X — ir/qp and tt/ (k F ± q F ). 

Unfortunately, since the intermediate temperature 
regime is absent for externally unscreened interactions, 
i.e., Tf « Tb, it appears to be rather difficult to decide 
which wavelength will eventually be more important for 
the experimental setup of Ref. |U| . For a well-screened 
interaction, based on the above discussion, we expect 
that the slow wavelength is more important. In that case, 
since Osdw-k can give a contribution to the ID spin den- 
sity s z {x) via the mechanism discussed in Sec. VIA , the 
s z correlations are dominated by the 2qp oscillatory part, 



(s z (x)s z (x')} ~ cos[2q F (x — x')] \x — x'\ 



h-2v 



For the distances 



where the correlations are not 



already vanishingly small, the cosine factor stays essen- 
tially at unity. Hence the SDW correlations show a pro- 
nounced ferromagnetic character. This reasoning offers 
an explanation for the ferromagnetic tendencies observed 
in Ref. jL8|. Our explanation could be experimentally 
checked by increasing the doping and hence q F . Even- 
tually, antiferromagnetism should be recovered for suffi- 
ciently strong doping. 



D. Superconductivity 

Predominant superconducting (SC) features have been 
predicted recently in two-chain models despite the repul- 
sive nature of the Coulomb interaction |4l],|^,|52| . For the 
nanotube, the inter-sublattice singlet and all triplet pair- 
ing operators cause irrelevant terms and therefore play 
no role. The dominant contribution then comes from the 
intra-sublattice singlet pairing (SCO) operator 
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OSCO ~ E P% aa %-a- 
paa 



This gives the bosonized form 



E 

rota 



~ cos[\/7r0c+] cos[\/7r0 c _] (6-14) 
x cos[v / 7T^ s +] sin[y / 7T'/'s-] — (cos <-* sin) . 

For T > 7],, this operator has scaling dimension 77 = (3 + 
l/K)/4 and is therefore irrelevant. In the intermediate 
temperature regime Tf < T < we find 77 = (3 + 
2/K)/8, and Osco is only subleading. However, for T < 
Tf, we obtain 



(O SCO (x)Osco(x')) ~\x-x'\ 



II-1/2K 



(6.15) 



Therefore singlet superconductivity becomes the domi- 
nant instability at extremely low temperatures and for 
externally screened interactions with K > 1/2. In prac- 
tice, superconducting correlations are thus of little im- 
portance in carbon nanotubes. 



E. Tunneling Density of States 

Under typical experimental conditions, the contact be- 
tween a SWNT and the attached Fermi-liquid leads is not 
always adiabatic. For nonadiabatic contacts, the physics 
of the conductance is related to single-electron tunnel- 
ing into the nanotube, which in turn is governed by the 
tunneling density of states (TDOS). A similar situation 
also arises in tunneling microscope experiments. Here the 
electron tunneling occurs close to the impurity position, 
the role of the impurity being played by the probe it- 
self. Therefore let us now discuss the TDOS of a SWNT. 
The presence of an open boundary, either due to the fi- 
nite length of the nanotube or due to a strong impurity, 
can alter correlation functions and, in particular, modify 
the anomalous exponents |35| ]. Hence, in calculating the 
TDOS close to the end of the SWNT, we have to employ 
the open boundary bosonization method. The latter can 
be developed for nanotubes in full analogy to spin chains 
p5| and quantum wires fl35| . Here we shall give the main 
ideas of the approach but avoid technical details, which 
the interested reader can find in Refs. |^5|,^5| . 

For simplicity, let us consider a semi-infinite nanotube, 
x > 0. The left-moving electrons in the nanotube then 
undergo perfect reflection at x = 0, where they are trans- 
formed into a right-mover going back to x — ► 00. Mathe- 
matically, the perfect reflection leads to boundary condi- 
tions imposed on the right- and left-moving electron field 
operators at x = 0, 



ip+aa{0) = exp(i(5 ) ip- a<T (Q) . 



(6.16) 



Here <5o is a scattering phase shift, which depends on 
the shape-confining potential. For example, if the end 
of the nanotube has dangling (but hydrogen-saturated) 
bonds, So is different compared to a closing cap, i.e., if 
the nanotube closes into a fullerene half-sphere. However, 
this phase shift is of no further interest here and can be 
chosen as zero by shifting the position of the boundary 



Due to the boundary condition (6.16), the right- and 
left-movers are not independent fields anymore. It there- 
fore makes sense to switch to a description in terms of 
only a right-moving field which is now defined for all x, 



*a<r(x) 



1p+aa(x) (X > 0) , 

ip- a(T (-x) (x < 0) 



(6.17) 



The advantage of this representation is that the right- 
moving field operator f M , being defined on the infinite 
interval, can be standardly bosonized, 



* QCT (x) 



1 



: exp[zyw(x)] , 



(6.18) 



where ip aa is a chiral Bose field. For brevity, we have 
omitted the Klein factors and the oscillating spatial expo- 
nent. Both cancel out in what follows. It is then helpful 
to define the chiral fields 



1 



4Vi 
1 

1 

1 



E 

eta 

E 

eta 

E 



^aa 1 



OOPa 



atp 



a a 



In terms of these fields, Eq. ( 6. IS ) takes the form 
* QCT = - 7 ^exp{i v / 7i : [*c+ + a<P c - +o-$ s+ +acr$ s _]} 

(6.19) 

If electron-electron interactions are neglected, the Qjs are 
free chiral fields with the Hamiltonian 



Ha 



^E r / dx ^La 9 x^ ra <r = 
rcca J ° 

/ dx & aa d x ^ aa = v / dx {d x $ jS ) 



(6.20) 



In the presence of interactions, however, the nonlocal rep- 



resentation (S.17) causes nonlocal interaction terms of 
the form p{x)p{— x) such that the electron densities are 
coupled at mirror-imaged points. Nevertheless, this in- 
teraction remains quadratic in the Bose field and can be 
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diagonalized. Following the steps of Ref. |3jJ , the appro- 
priate Bogoliubov rotation is 



<I>c+(x) -► C^ c+ (x) - S$ c+ (- 



(6.21) 



where c = cosh </>o and s = sinh 0o parametrize the Lut- 
tinger liquid exponent K — exp(20o)- 

Two remarks are in order here. (1) We only consider 
the Luttinger liquid phase T > Tf, in the following. Then 
the only important interaction term is H^p S . In the low- 
temperature massive phases, the single-electron TDOS 
develops a gap, but there could still be gapless higher- 
order multi-electron scattering processes. (2) We con- 
sider the zero-momentum limit for the charge exponent 
A, i.e., distances (times) larger than the Coulomb screen- 
ing length. The behavior of the correlation functions 
for distances closer to the boundary than the screening 
length is fairly complex, see Appendix C of Ref. [jj5| . 

Combining Eqs. (6.19) and (6.21), we obtain the main 
formula of the open boundary bosonization, 



: exp {iy/ir [c$ c+ (x) - s$ c+ (- 



+a$ c _(a;) + ct$ s+ (x) + aa$ s - (x)]} 



(6.22) 



It follows that close to the boundary, i.e., for max(x, y) <C 
vt, the single-electron Greens function is of the form 



(&(x,t)V(y,0)) „ t -(i/*+3)/4 



(6.23) 



The boundary scaling dimension of the electron field op- 
erator is therefore 



* 1 3 
A= 8A + 8' 

as opposed to its bulk scaling dimension, 



A = i 

16 \K 



(6.24) 



(6.25) 



Making use of the text-book definition of the density of 
states as the imaginary part of the single-electron Greens 
function, we obtain for the TDOS at the end of the nan- 
otube, 



2A-1 



(6.26) 



which is different from the bulk density of states 

PbuikM ~ w 2 *- 1 . (6.27) 

Since A > 1/2 for A' < 1, the TDOS always vanishes 
with the frequency (or the temperature) approaching 
zero, similar to the situation encountered in quantum 
wires [fj3| . For a Fermi liquid, we would instead have 
A = A = 1/2, and the exponents in Eqs. (3.26) and 



(S.27) are both zero. For the (10, 10) nanotube of Ref. 



fl8|| , taking our previous estimate A ~ 0.18, the expo- 
nents in Eqs. ( 6.26| ) and (3.27) are 1.13 and 0.46, re- 
spectively. Therefore single-electron tunneling is strongly 
suppressed close to the end of the nanotube. Finally, we 
stress again that these results are also valid in the vicin- 
ity of a strong impurity. At low energy scales, the SWNT 
is effectively cut into two independent parts because of 



the impurity, see Sec. VII 



VII. TRANSPORT AND CONDUCTANCE LAWS 

In this section, we discuss transport through a SWNT. 
Since non-Fermi liquid laws are only pronounced if con- 
tact resistances between the attached leads and the nan- 
otube are smaller than h/e 2 , we assume good con- 
tact to the transport voltage sources. Since there are 
four transport channels (j5), a SWNT adiabatically con- 
nected to external leads exhibits the perfect conductance 
G = Go = 4e 2 //i in the absence of impurity backscat- 
tering. Even in the presence of weak backscattering, G 
approaches Go at sufficiently high temperatures. 

The effect of impurities on the conductance can be very 
pronounced in a Luttinger liquid J33| . In a SWNT, there 
are several possible sources for impurities as discussed in 
the Introduction. The dominant coupling to the impurity 
is due to the CDW7T operators 0\ and 02 in Eqs. (6.4) 
and (6.5). Since both operators have essentially the same 
effect, we focus on 0\ and get the form 



iAmp = J dx m(x){cos[v / 7r#c+ + ^qpx] cos[t/tt9 c -] 

x sin[\/7r0 s+ ] sinfv^s-] + (cos <-> sin)} . (7.1) 

From the perturbation series in m{x) for the partition 
sum, one can check that the higher-order operators of 
Sec. VI B| are implicitly contained in Eq. (7.1). There is 
no need to explicitly take them into account in H lmp . 

Let us analyze the consequences of this Hamiltonian for 
the simplest case of a single pointlike impurity, m(x) — 
mS(x). We assume an impurity without any internal dy- 
namics, e.g., a frozen twist of the nanotube. The con- 
ductance corrections SG defined by G = Go — SG depend 
on the temperature regime under study. Perturbation 
theory yields in order m 2 the power laws 



5G ~ T {K - l)/2 (T > T b ) , (7.2) 
~ T (2K ~ 5 ^ 4 (T f <T< T b ) . (7.3) 

At very low temperatures, T < Tf, one is in the strong- 
coupling regime and perturbation theory does not ap- 
ply anymore. In fact, we find that the conductance is 
totally suppressed at zero temperature. In the re gime 
T < Tf, the average value of the impurity operator (7.1) 
vanishes, while its correlation function decays exponen- 
tially in time. At first sight, the impurity has therefore 
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no effect. However, it may still generate a relevant op- 
erator. In order to see how this happens, we calculate 
the second-order correction to the action due to the local 
impurity operator, 

SS 2 ~m 2 J dr J dr' exp[- 7 |T - r'|] 

x cos[V7r0 c+ (O,r)] cos[-y/n9 c+ (0,r')] , 

where 7 ~ rrif characterizes the decay rate of correla- 
tions of the gapped degrees of freedom involved in the 
impurity operator (|7.l|). Since we are interested in the 
low-energy (long-time) behavior of the system, we are 
eligible to contract the r and t' variables in the above 
integration. Then, using the operator identity ( |5.9| ), the 
total charge field satisfies the effective Hamiltonian 

H cS [6 c+ ] = H [6 C+ ] + A cos [Virf 9c+ (o)] , (7.4) 
ira 

where the effective impurity strength is estimated to be 



A 



~/rrif, and we have rescaled 0, 



c+ 



K0, 



<:+■ 



Of 



course, the exponent found in Eq. (|7.5| ) below is indepen- 
dent of A. Note that the result ( [7.4[ ) also follows from 
Eq. (6.7) for T < Tf by virtue of the pinning conditions 
discussed in Sec. 



The effective Hamiltonian (|7.4| ) is now identical to the 
one describing an impurity in a spinless Luttinger liquid. 
The reasoning of Kane and Fisher |}3) therefore fully 
applies and yields the low-temperature conductance 



(7.5) 



This implies a total suppression of transport through the 
nanotube at T = in the presence of a single, arbitrarily 
weak impurity. 

Finally, we mention that Luttinger liquid behavior 
could also be observed in experiments on crossed nan- 
otubes. Noncquilibrium transport through crossed Lut- 
tinger liquids has been studied theoretically in Ref. |36| . 
Following this analysis, a distinct nonlinear dependence 
of the current through one SWNT (j — 1) on the cross 
voltage applied to the other SWNT (j = 2) can be ex- 
pected. Assuming a point-like contact between both nan- 
otubes at x = 0, the only important coupling mechanism 
is of electrostatic nature, 



H c = Ap x (0)p 2 (0) 



(7.6) 



which becomes relevant if the scaling dimension r\ of 
the individual local density operators Pj=i,2(0) satisfies 
7] < 1/2. This condition can only be met by the higher- 
order CDW7T contributions in Sec. 



VI B 



In the Luttinger 

liquid regime T > Tf,, the leading operator is 0\, DWlx in 
Eq. ( 6.10| ), with i] — AK. The coupling (7.6) is then irrel- 
evant unless the condition K < 1/8 is fulfilled. According 



to Eq. (3.23), such strong correlations could be achieved 
by studying sufficiently long nanotubes. Furthermore, 



at lower temperatures T < T;,, the dominant coupling 
comes from 0\ in Eq. (6.7), with 77 = K. Therefore the 



relevancy condition is now less stringent, K < 1/2, and 
the effects predicted in Ref. [|6| should show up already 
for short nanotubes at sufficiently low temperatures. 



VIII. IMPURITY SCREENING 

In this section, we study the electronic screening cloud 
induced by an impurity sitting at, say, x = 0. This 
could be a substitutional atom replacing one of the car- 
bon atoms on the tube surface (like B or N), an interca- 
lating atom (like Na) sitting inside the nanotube at some 
radius r < R, or an incident positron. The (externally 
unscreened) Fourier transformed potential for impurity 
charge Ze is [B8| 



V imp (q) = 4irZe 2 I (qr)K (qR) , 



(8.1) 



with the modified Bessel functions Iq and K$ j46|. The 
small-g response of the 7r electrons to this perturbation 
can be computed in an exact manner by linear screening, 
and the RPA treatment and the corresponding results of 
Ref. fully apply. However, for the finite wavevectors 
listed in Eq. ([O]), linear screening is known to break 
down in one dimension The ID Friedel oscilla- 

tion cannot be computed by an RPA-like treatment. 

The impurity strength m determining the amplitude 
of the Friedel oscillation is given by V r i m p(2gp') [for wave- 
length A = ir/qp], or Vi mp (2kp) [for A = n/kp — 
■k /{kp ±qp)]. In the latter case, m is strongly dependent 
on the position of the impurity. If the impurity sits in the 
center of the nanotube (r = 0), we get exponential sup- 
pression of m due to Eq. ( |S.lD . This implies a strongly 
reduced amplitude of the Friedel oscillation. However, 
the actual value of m does not affect the power laws re- 
ported below. These are universal, i.e., independent of 
strength, position r, o r nat ure of the impurity. 

As discussed in Sec. VI A , the ID density operator q(x) 
has several contributions. The expectation value of the 
slow part p(x) in the presence of the impurity can be 
obtained by RPA j38| and is not further discussed here. 
The contributions beyond p(x) are due to the CDW or- 
der operators in Sec. VL These arise from mixing the 
different types of ID fermion operators. 

Following Ref. [Q , it is straightforward to extract the 
power law decay of the Friedel oscillation. For T >T\, 
the CDW-7T operators Ox and O2 specified in Eqs. (6.4) 
and (6.5) lead to the simultaneous presence of the wave- 
lengths ir/qF and 7r/(fcp ± qp). Contrary to the RPA 
analysis [p8[ , our exact treatment of Coulomb interac- 
tions yields for x ^> a • 



(q2q F (x)) ~ COS 

~ cos 



[2q F x] (x/x y {3+K)/4 (x » 10) 
[2q F x] (x/x )- (1+K)/2 (x « x ) 
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(q2k F (z)) ~ cos[2(fc F ± q F )x] {x/x )- {S+K)/4 (a; » x ) 
~ cos[2(fc F ± q F )x] (x/x y (1+K)/2 (x < x ) . 

We omit possible phase shifts here. Furthermore, xq sets 
the appropriate length scale, which can be different for 
the various oscillation periods since it depends on the 



impurity strength, xq ~ m 



-4/(l-K) 



The asymptotic 



exponents for i>io can be ob tained from open bound- 
ary bosonization, see Sec. VIE, while the behavior close 



to the impurity follows from perturbation theory in m. 
Here the important point is that the Friedel oscillation is 
always slower than the standard 1/x Fermi liquid result. 
Furthermore, the decay becomes slower as the impurity 
is approached. Following the reasoning of Ref. fjql , the 
slow decay of the Friedel oscillation is the physical rea- 
son for the conductance suppression by a single impu- 
rity discussed in Sec. VII. Interestingly, we have again 



several different oscillation periods, namely a slow one, 
A = Tr/qF, which is superimposed onto rapid oscillations 
with wavelength A = ir/(kp ± qp). The Friedel oscil- 
lations induced by a strong impurity or by the ends of 
the nanotube can be quite pronounced and should be 
detectable by placing a STM tip close to the nanotube. 
The power laws (8.2) and (|S.3D are modified due to 



the 8qF operator (6.1C) in the case of extremely strong 
correlations. This operator plays no role at temperatures 
T < Tb, but in the Luttinger liquid regime it becomes 
dominant as soon as K < 1/5. Instead of Eqs. (S.2) and 
(3.3), we then find the wavelength ir/Aqp, 



(q8q F ( x )) ~ coa[8q F x] (x/x ) 4K (x^>x ) 
~ cos[8qFx] | ln(a;/xo)| (x <C Xq) 



(8.4) 



The only logarithmically slow decay of the Friedel oscil- 
lation close to the impurity follows from Ref. pj] ]. 

The behavior of the Friedel oscillations discussed so 
far applies only to the "Luttinger liquid" regime, T > 
Tjj. In the intermediate temperature regime, Tf < T < 
Tb, and for not too strong correlations, 1/2 < K < 1, 
the exponents in these power laws become — (3 + 2K)/8 
for x ^> xo, and (1 — 2K)/4 for i < i , respectively. 
For strong correlations, K < 1/2, the squared CDW7T 
operators of Sec. VI B are more important than 0\ and 
02. In that case, only the slow oscillation period ir/2qp 



corresponding to Eq. (6.7) shows up, and we obtain 



(qAq F (x)} ~ cos[4q F x] (x/xq) K (2; > xo) 
~ cos[4g F a;] | ln(x/xo)\ (x -C Xq) 



(8.5) 



Finally, at the lowest temperatures, T < Tf, the same re- 
sult is found for K < 1/2. For K > 1/2, the perturbative 



exponent 



„1-2K 



applies to the short-distance behav- 



.5) 



ior, with the same asymptotic exponent as in Eq. 

The overall behavior of Friedel oscillations in nan- 
otubes is considerably more complicated than predicted 
by linear screening. We find several different slow alge- 
braic power laws in various regions of parameter space. 



.Slj'hey are all slower than the standard 1 jx decay. Further- 
more, there exist different wavelengths for the respective 
dominant contributions. If measurements of the charge 
screening cloud become feasible, one could extract the 
important correlation parameter K and directly observe 
the various regions of the phase diagram corresponding 
to Table Q from the Friedel oscillation. 

IX. CONCLUSIONS 

In this work we have formulated and analyzed the ef- 
fective low-energy theory of isolated single-wall carbon 
nanotubes. Long and clean nanotubes have already been 
fabricated, and we believe that the non-Fermi liquid ef- 
fects discussed here will be observed in the near future. 

In fact, there is some evidence that one of these ef- 
fects, namely the predominance of ferromagnetic correla- 
tions, has already been observed in Ref. |fl8|| . Other sig- 
natures of non-Fermi-liquid behavior include anomalous 
interaction-dependent power laws for the conductance or 
the tunneling density of states. Nanotubes are promis- 
ing candidates for revealing the Luttinger liquid behavior 
generally expected for one-dimensional metals. Accord- 
ing to our theory, at temperatures above the scale Tb 
(where Tb ~ 0.1 mK for the setup of Ref. a nan- 

otube behaves as a spin-i Luttinger liquid with an addi- 
tional flavor index. At lower temperatures, more compli- 
cated phases characterized by gaps in the neutral modes 
emerge. 

A special feature of doped nanotubes is the simultane- 
ous occurrence of different wavelengths for the various or- 
der parameter correlations. The doping is characterized 
by a small wavevector qp <C hp, and one quite generally 
finds pronounced instabilities for the long wavelengths 
A = ir /qp,Tr /2qp and ir/Aqp, besides the short wave- 
length A = 7r /{hp ± qp). Another interesting aspect of 
nanotubes is the length (L) dependence of their proper- 
ties. For instance, considering a SWNT with sufficiently 
small L, our continuum approach for explaining ferro- 
magnetic tendencies should eventually be replaced by a 
Hund's rule applying to the "molecule" . In our opinion, 
this transition from a molecule to the "solid-state" limit 
of large L deserves further attention. 

The physics of nanotubes is likely to reveal further sur- 
prises in the future. Besides the exciting prospect of find- 
ing Luttinger liquid behavior in a potentially very clean 
way, nanotubes might act as basic elements for molecu- 
lar electronics devices, e.g., as highly conducting wires. 
Another line of development could employ networks of 
nanotubes, where novel correlation effects can be antici- 
pated. We hope that our theory will be useful in address- 
ing and resolving these issues. 
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TABLE I. Dominant correlations in an armchair SWNT 
as a function of temperature and of the correlation parameter 
K < 1, with the respective scaling dimension r\ and wave- 
length A. 
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FIG. 1. Schematic bandstructure of a metallic SWNT. 
A right- and left-moving branch (r = ±) is found near each 
of the two Fermi points a = ± corresponding to K and K' , 
respectively. Right- and left-movers arise as linear combina- 
tions of the sublattices p = ±. The Fermi energy (dashed 
line) can be tuned by an external gate. 



(a) aFS 



p aa p aa 




\>'ol'o' p'a' 'a' 



(b) aBS 



p aa p -aa 




\>'-aa' \>'ao' 

FIG. 2. Allowed interaction processes away from 

half-filling, (a) Forward scattering (aFS). The incoming ID 
fermions do not change their quantum numbers p, a, a dur- 
ing the scattering event, (b) Backward scattering (aBS). One 
ID fermion is scattered from Fermi point a to the opposite 
Fermi point — a, with the other being scattered from —a to 
a. Straight lines denote ID fermion propagators, dashed lines 
the effective ID Coulomb interaction. 




FIG. 3. Microscopic arrangement of carbon atoms around 
the waist of an armchair SWNT (shown for n = 4). Circles 
and crosses denote the two distinct sublattices p = ±, and 
d = a/V3. 
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FIG. 4. RG flow of the coupling constants for / = 0.025 



and b = .05. Units are such that a = 1. Following Eqs. (2.25) 
and ( |2.29| ), these initial values apply to a (10, 10) SWNT with 
externally unscreened Coulomb interactions. 
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FIG. 5. Same as Fig. W for / = 0.05 and b = 0.025. 



23 



